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Abstract 

In earlier papers, the quantum amplitude for processes involving the formation and 
evaporation of black holes was calculated by means of a complex-time approach. Instead 
of taking a more familiar approach to black-hole evaporation, we simply followed Feyn- 
man's + is approach in quantum field theory. The Lorentzian time-interval T , measured 
at spatial infinity between a pair of asymptotically-flat space-like hypersurfaces Ej and E_p 
carrying initial and final boundary data for the gravitational and other fields, is rotated: 
T — ► \T\ exp(— iS) , where < 5 < n/2. Classically and quantum-mechanically, this 
procedure is expected to lead to a well-posed boundary-value problem. Thus, what we 
have done is to find quantum amplitudes (not just probability densities) relating to a pure 
state at late times following gravitational collapse of matter to a black hole. Such pure 
states, arising from gravitational collapse, are then shown to admit a description in terms 
of coherent and squeezed states. Indeed, this description is not so different from that aris- 
ing in a well-known context, namely, the highly-squeezed final state of the relic radiation 
background in inflationary cosmology. For definiteness, we study the simplest model of 
collapse, based on Einstein gravity with a massless scalar field. Following the complex 
rotation above, one finds that, in an adiabatic approximation, the resulting quantum am- 
plitude may be expressed in terms of generalised coherent states of the harmonic oscillator. 
A physical interpretation is given; further, a squeezed-state representation follows. 

PACS numbers: 04.70.Dy, 04.60.-m, 03.65 

1. Introduction 

The treatment given here of black-hole evaporation depends on Feynman's + ie approach 
[1]. We begin by reviewing this approach, which in [2-8] was described and applied to 
quantum amplitudes (not just probabilities) for particle production, following gravitational 
collapse to a black hole. In specifying either the classical boundary-value problem or the 
quantum amplitude to be computed, one takes an initial space-like ypersurface Ej and a 
final space-like hypersurface E p . For simplicity, imagine that we are in the asymptotically- 
flat context and, further, that each of E/ and E p has the topology of IR 3 . One then 
specifies suitable asymptotically-flat boundary data for Einstein gravity and for any other 
fields present on Ej and E_p . If there are no fermionic fields present in the Lagrangian, 
then one might expect the gravitational boundary data to consist of the intrinsic positive- 
definite (Riemannian) 3-geometry, given by the intrinsic 3-metric hij = gij = 1,2,3) 
on a surface {x° = const.}. Here, g^ v (/i, v = 0,1,2,3) gives the 4-dimensional 'space- 
time' metric. For a complete specification of the problem (if indeed such a boundary-value 
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problem is well-posed - see below), one further needs to give the Lorentzian proper- 
time interval T which separates Ej from £p , as measured at spatial infinity. In the 
papers [3,4], as in the present paper, we take the simplest possibility for the matter fields 
present, namely, a real massless scalar field cf> . Then, suitable boundary data on Ej and 
on Ep are expected to be For further simplicity, without great conceptual 

loss, we assume here that the initial data on Ej are spherically symmetric, corresponding 
to a diffuse slowly-moving initial matter distribution, prior to gravitational collapse. The 
final gravitational and scalar data are taken to have a 'background' spherically-symmetric 
part, together with small non-spherical perturbations, which correspond to gravitons and 
massless-scalar particles; such data represent a possible field configuration on a space-like 
slice of topology M 3 at late times, for which a quantum amplitude should be calculable. 

In Feynman's +ie approach [1], one rotates the time-interval T into the lower-half 
complex plane: T — > |T| exp(— id), where < 5 < tt/2. Except in the extreme case 
5 = tt/2, one has classically a complex boundary- value problem, with a complex classical 
4- metric and scalar field <fi m the interior, for the given (real) boundary data (hij , <f>) on 
Ej and E^ . In the special case 5 = tt/2 , one has a real Riemannian time separation \T\ 
at spatial infinity, and therefore one expects to have a well-behaved real elliptic boundary- 
value problem (modulo gauge) for (g^ , <p) . From comparison with the case of linearised 
massless scalar fields [4,9], one expects that the complex case (0 < 5 < tt/2) should be 
strongly elliptic [10]; that is, despite the complex boundary parameter T, the boundary- 
value problem should continue to be well-posed, with analytic solutions, existence and 
uniqueness. By contrast, the 'boundary- value problem' for 5 = , in which one asks for a 
real Lorentzian-signature solution to the boundary-value problem, is well known to be ill- 
posed [11]: it is inappropriate to pose boundary- value problems for hyperbolic (wave-like) 
systems. 

For the black-hole collapse problem, we consider boundary data of the type described 
above, with weak-field linear-order anisotropic perturbations (A , A^ 1 )) away from 
spherical symmetry on Ep, but with no anisotropy on £ j , where A is a small parameter. 
The Lorentzian classical action of the exactly spherically-symmetric solution (here assum- 
ing 5 > ) described by g$ = 7 MZ , and 4>^ = $ , will be written as S^iass • Given the 
small (but non-zero) anisotropy of the final data and hence of the interior classical solution 
(g^v > 4>) ? one considers the asymptotic expansions: 

9ixv ~ lixv + \g$ + \ 2 g$ + ... , 
(j) ~ $ + A(/> (1) + A 2 (/> (2) + ... . 

For such linearised boundary data, with corresponding perturbative expansions (1.1) for 
the classical solution, the classical Lorentzian action Sdass will have an asymptotic expan- 
sion of the form: 

Sclass ~ ^class + ^cliss + ^class + • • • • (I- 2 ) 

Here, the perturbative factors A 2 multiplying S^l ss , A 3 multiplying S^fl ss , etc., have been 
omitted, in order to simplify the notation later. The first correction S^ ag is the second- 
variation classical action and is bilinear in the linear-order [Xh^ , X(p^) corrections 
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to the boundary data on T, F . As in [7], one can indeed evaluate S^iass j°^ nt ^ as a 
functional of the linearised boundary data (A , A 4>^) F an d a function of the complex 
variable T . One then computes the corresponding semi-classical quantum amplitude, 
proportional to exp^S^j^s) , and one can also include loop corrections, if appropriate. In 
fact, it appears likely, first, that any quantum field theory which includes gravity must be 
invariant under local supersymmetry (and hence a theory involving supergravity) if it is 
to have meaningful quantum amplitudes, and, second, that such a quantum theory may 
actually have finite quantum amplitudes, even though these may include loop corrections, 
unless (for example) the theory is pure N = 1 supergravity without supermatter [9,12]. 
In that case, the quantum amplitudes considered in the black-hole problem can certainly 
be truncated as exp^S^ 2 ^), provided that the frequencies involved in the boundary data 

(A/i^j , A^*- 1 -*) are well below the Planck frequency. Finally, then, following Feynman's 
+ is prescription in the present context of black-hole quantum evaporation, one recovers the 
Lorentzian quantum amplitude (again, not just the probability density) for the quantum 
state including (say) created particles present at late times, by taking the limit of the 

semi-classical amplitude exp(zS'^ ss ) as 5 — > 0+ . 

As seen in [4,6,8], the black-hole radiation has the usual thermal spectrum, at the 
temperature 1/8% Mj . Unless otherwise stated, we employ Planckian units, taking: ks = 
c = h = G = 1 . Here, Mj is the ADM (Arnowitt-Deser-Misner) mass of the 'space- 
time' [27]. As described in [3] and References therein, the ADM mass on the initial 
surface Ej must equal the ADM mass on the final surface Ep , in order that the classical 
boundary- value problem should be well-posed, and in particular that the space-time metric 
components should have the expected fall-off properties at spatial infinity. In the slightly- 
complexified regime, where the (large) time-interval T at spatial infinity obeys Im(T) < 
, the geometry is accurately approximated by a Vaidya metric [6,20], with slowly- varying 
mass function m{u) , where u is a 'retarded-time' coordinate. If one wished, one could 
pursue the (classical) perturbation theory of [6] further, so as to describe accurately the 
slightly-complex classical geometry, given anisotropic weak-field final data on the late-time 
hypersurface T,p , whose mass agrees, as above, with the mass of the initial-data set on 
Ej • That is, not only would one arrive at the usual retarded-time dependence of m{u) 
for a radiating black hole - one would also (given sufficient labour) calculate the detailed 
evolution backwards in time of the weak-field anisotropic perturbations. 

Of course, the usual definition of a black hole depends on the space-time metric being 
real, of Lorentzian signature. In contrast, the classical Einstein boundary- value problem 
is only expected to be well-behaved when Im(T) < , for a time-interval T measured at 
infinity. Thus, strictly speaking, it is inappropriate to use the term 'black hole' in relation 
to the complexified classical solution of the previous paragraph. However, it would be a fair 
use of language to say that the infilling Vaidya-like classical solution for the boundary- value 
problem of the previous paragraph is a 'black-hole intermediate state'. 

Readers who are accustomed to the original Lorentzian-signature derivation of black- 
hole evaporation will be used to the notion of radiation (scalar- field, Mazwell, etc.) piling 
up around the future event horizon, and then undergoing an enormous redshift, depending 
in a specific way on the mass Mj , as the radiative fields move out towards future null 
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infinity. The detailed form of this redshift is intimately connected with the thermality and 
temperature of the black-hole radiation. But this radiative behaviour can be learnt equally 
well through study of high-frequency or WKB solutions of the separated wave equation, 
as in Eqs. (2.4,5) below. The WKB transmission and reflection coefficients describe the 
thermal radiation. And WKB investigation of wave equations such as Eqs. (2.4,5) is also at 
the base of the present approach [4,6,8]. One might say that detailed knowledge concerning 
the (Lorentzian) event horizon is 'imbedded' in the relevant spin-s wave equations, such as 
Eqs. (2.4,5). Thus, it should not be surprising that one can arrive at the thermal spectrum 
using either geometrical considerations of the Lorentzian horizon or analytic arguments 
requiring a complex metric (for which case the event horizon is undefined). In similar 
fashion, one may still speak here, if desired, about particles which fall into the hole as well 
as particles which travel out to infinity; these correspond to the usual basis of two WKB 
solutions of the radial wave equation. 

In the present paper, we study quantum amplitudes found via Feynman's approach, as 
discussed and calculated in [2-8], but now in the context of coherent states [13], which re- 
semble 'classical states', and of squeezed states [14], which are purely quantum-mechanical. 
As above, our motivation originated with the study of the final radiation which remains 
after a black hole has evaporated completely, but there are strong connections also with the 
relic Cosmic Microwave Background Radiation (CMBR) induced by inflationary cosmo- 
logical perturbations. In fact, particle creation by black holes has several similarities with 
cosmological particle creation, despite the lack of asymptotic flatness in the cosmological 
case. Cosmological and black-hole particle creation both require a time-dependence in the 
metric. This in turn has led to descriptions, in terms of coherent and squeezed states, of 
quantum phenomena in curved space-time; some of the earlier examples include [15-17]. 

In inflationary cosmology, the field modes are in their adiabatic ground state, with 
short wavelengths near the start of inflation. This is related to the assumption that the 
universe was in a maximally-symmetric state at some time in the past [18,19]. Due to the 
accelerated expansion of the universe during inflation, quantum fluctuations are amplified 
into macroscopic or classical perturbations. The early-time fluctuations lead to the for- 
mation of large-scale structure in the universe, and also contribute to the anisotropies in 
the CMBR. The final state for the perturbations is a two-mode highly- squeezed state for 
modes whose radii are much greater than the Hubble radius, with pairs of field quanta 
(having opposite momenta) being produced at late times [16,17]. Tensor (s = 2) fluctua- 
tions in the metric, for example, are predicted to give rise to relic gravitational waves. On 
the other hand, electromagnetic waves (s = 1) cannot be squeezed in the same way during 
the cosmological expansion, because they do not interact with the external gravitational 
field in the same way. 

In either the cosmological or the black-hole case, one works within an adiabatic ap- 
proximation for the perturbative modes. Writing k for a typical perturbative frequency, 
one requires k ^> aH in the cosmological case, where H = (a/a) , with a(t) denoting 
the scale factor. In the black-hole case, the space-time geometry at late times, in the 
region containing a stream of outgoing radiation, is given by a Vaidya metric [20-21] with 
a slowly- varying 'mass function' m(t, r) . The adiabatic condition then reads k 3> \m\/m. 
Indeed, for an evaporating black hole, at all except the last moments of evaporation, the 
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frequencies of interest in the evaporating modes do typically exceed \m\/m, namely, the 
inverse time-variation scale for the black-hole mass. In other words, the period of a wave 
of interest is typically much smaller than the time-scale for variations in the background 
gravitational field. The black hole interacts negligibly with the emitted particles, and the 
time between successive emissions is comparable with the black-hole mass [22] . 

In the cosmological physical description of the corresponding phenomenon is 

that, when the wavelength is comparable with or larger than the Hubble radius, ampli- 
fication of the zero-point quantum fluctuations takes place. As a further aspect of this 
analogy, the redshifting of the radiation in the black-hole background space-time is deter- 
mined by the total mass Mj , and correspondingly by the Hubble parameter H~ x in the 
cosmological case. 

In this paper, we apply the squeezed-state formalism to black-hole evaporation, while 
maintaining a comparison with inflationary cosmology - see [23] for further comparisons. 
In the case of inflationary cosmology, the quantum evolution of cosmological perturbations 
(density, rotational and gravitational), which begin in an initial vacuum state, follows 
essentially a set of Schrodinger equations [19]. The state of the perturbations is transformed 
into a highly-squeezed vacuum state, with many particles, having a large variance in their 
amplitude (particle number), but small (squeezed) phase variations. The squeezing of 
cosmological perturbations may be suppressed at small wavelengths, but it should be 
present at long wavelengths, especially for gravitational waves [24]. These perturbations 
also induce the anisotropies at large angular scales, as observed in the CMBR. Their 
wavelengths today are comparable with or greater than the Hubble radius. The above 
amplification of the initial zero-point fluctuations gives rise to standing waves with a fixed 
phase, rather than travelling waves. The relic perturbations in the high-squeezing or WKB 
limit can be described as a stochastic collection of standing waves [16,17]. This paragraph 
has reviewed the cosmological will be seen below, a similar picture emerges in the 

application to black-hole evaporation. 

In Sec. 2, we outline the main features of the above complex approach to the calcula- 
tion of quantum amplitudes (not just probabilities) for perturbative data (spins s = 0, 1, 2) 
prescribed on a late-time final hypersurface E p . For this procedure to be well-posed, one 
has first to rotate: T — > \T\ exp(— id) into the lower half-plane. The resulting stand- 
ing waves, which originate from setting Dirichlet boundary data on the initial and final 
space-like hypersurfaces, turn out to correspond to a highly-squeezed final state for late- 
time black-hole radiation. In the adiabatic approximation, the fixed phases correspond to 
discrete frequencies in the remnant (quantum) radiation from the evanescent black hole. 

Sees. 3, 4, 5 describe coherent states, generalised coherent states and squeezed states, 
respectively. In Sec.6, the small angle 5 , through which the time T at spatial infinity is 
rotated into the complex: T — > |T| exp(— iS) , is related to the large amount of squeezing 
which has been applied to give the final state. We also discuss briefly the normalisation of 
the probability density, and demonstrate the existence of large peaks and troughs in the 
spectrum of the radiation reaching future null infinity, due to the standing-wave pattern 
of the perturbations. A short discussion of entropy and squeezing is given in Sec. 7, and 
possible classical predictions are considered in Sec. 8. Sec.9 contains a brief Conclusion. 
Some technical results are given in the Appendix. A briefer account of this work has 
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appeared in [25]. 

2. The quantum amplitude for late-time data 

Consider first the effect on the classical boundary-value problem of a rotation into the 
complex of the time-interval T by a moderately small angle 5 , as above. The resulting 
classical solution [g^ v , <j>) of the coupled Einstein/massless-scalar field equations will be 
somewhat complexified, by comparison with a Lorentzian-signature solution. By suitable 
choice of coordinates (£, r, 9, <p) , the spherically-symmetric 'background' part of the metric 
may be written in the form [2,4] 

ds 2 = - e b dt 2 + e a dr 2 + r 2 (d9 2 + sin 2 6> dip 2 ) , (2.1) 

where b = b(t , r) , a = a(t , r) , and the spherically-symmetric 'background' part $ of 
the scalar field has the form $ = <£>(£, r) . The spherically-symmetric functions a, b and 
$ must, of course, be complex- valued. The coupled Lorentzian-signature Einstein/scalar 
field equations for this spherically-symmetric configuration are given by the analytic con- 
tinuation of the Riemannian field equations Eqs. (1.9-13) of [7], on making the replacement 

t = rexp(-n?) , (2.2) 

where r is the 'Riemannian time-coordinate' of [7], and where the real number $ should 
be rotated precisely from to tt/2 . 

The small anisotropic perturbations in the boundary data on the final late-time hyper- 
surface T,p consist, in the language of Sec.l, of the perturbed part A/ijj'' of the intrinsic 
3-dimensional spatial metric hijF on Sf , together with the perturbation A^ 1 ) of the 
scalar field 4>f on Hp- The classical solutions resulting from these perturbed boundary 
data correspond to gravitons and to massless-scalar particles, propagating on the (com- 
plex) spherically-symmetric classical background [g^ , $) . For example, the field in 
the interior of the space-time may be decomposed as in Eq.(6) of [3]: 

^ oo m=£ 

(f>^(t,r,e t( p) = -J2 Yl ^-(°) Rim(t,r) . (2.3) 

£=Q m=-£ 

Here, Ye m (Q) denotes the (£, m) scalar spherical harmonic of [23]. The scalar field equation 
decouples for each (£ , m) , leading to the mode equation 

(e( & - a )/ 2 d r ) R tm - {d t ) 2 R im - - (d t {a-b)) [d t R lm ) - V e (t,r)R £m = , (2.4) 
where 

e 6( *' r) ( tt(ll 2m(t,r) 



V t {t,r) = —j-^+IJ-h—L^J (2.5) 

would be real and positive in the Lorentzian-signature case. The 'mass function' m(t, r) , 
which would equal the constant mass M/ for an exact Schwarzschild geometry [27], is 
defined by 

e -a(t,r) = 1 _ (26) 
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A corresponding angular decomposition can be given for weak gravitational-wave pertur- 
bations about the spherical background [7,28-30]. 

In most regions of the classical space-time, except for the central region where the 
black hole is formed, the metric functions a(t, r) and b(t, r) must vary only 'slowly' or 
'adiabatically'. This allows one to study, in such a region, a radial mode solution for (say) 
a perturbed scalar field, in which a further separation of the time dependence can be made 
approximately [2,25]: 

Re m (t,r) ~ exp(ikt) £ktm(t,r) . (2.7) 

Here, £k£m(t, r) varies 'slowly' with respect to t . This will occur near spatial infinity and 
also, provided that the time-interval T is sufficiently large, in a neighbourhood of the final 
hypersurface T,p . The mode equation (2.4,5) then reduces [4] to 

e^ 2 £ (e«- a V* + (k 2 - V ( ) Hum = . (2.8) 

As seen in [2,6], the spherically-symmetric background metric in this region can be rep- 
resented to high accuracy by a Vaidya metric, which describes the outflow of massless 
matter, which is spherically symmetric, on the average. The principal condition for the 
validity of the adiabatic expansion is [6] that 

\k\ > ^ . (2.9) 
m 

In studying the behaviour of solutions of the radial mode equation (2.8), it is natural 
to define a generalisation r* of the standard Regge- Wheeler or 'tortoise' coordinate r* s in 
the Schwarzschild geometry [27], by 

° = e^/ 2 ^- . (2.10) 



dr* dr 
In the exact Schwarzschild metric, this gives 



r 

s 



r + 2Mj log^(r/2M/)-l^ . (2.11) 



The approximate (adiabatic) mode equation (2.8) then reads 



dr 



k€m 



*2 



(k 2 - V^j i Mm = . (2.12) 



The procedure adopted in [2,6] involves choosing a convenient set of suitable radial 
functions {iktm{r)} on the final surface Tip , since it is here that the non-trivial boundary 
data are posed. The mode equation (2.12) does not depend on the quantum number m, 
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whence we may choose £k£m( r ) = Cki( r ) , independently of m. The boundary condition 
of regularity at the spatial origin {r = 0} [4,6] implies that 



£ u ( r ) = constant x (kr) i+1 + o((£r) £+3 ) (2.13) 

as r — > + . To examine the boundary condition on the ikiix) as r — > oo , note that the 
potential Vt(r) decreases sufficiently rapidly, as r — > oo , that a real solution to Eq.(2.12) 
behaves near {r = oo} according to 

€ki(r) ~ (z M exp(ikr* s ) + z* M exp(-ikr* s )^j . (2.14) 

Here, the Zki are certain dimensionless complex coefficients, which must be determined via 
the differential equation (2.12) together with the regularity conditions. Further, there is a 
natural normalisation of the basis {£ki( r )} of radial wave-functions, as discussed in detail 
in [6]. 

Given an appropriately normalised basis {Cki(f)} of radial wave- functions on the 
final hypersurface S p , one can expand out the interior linearised classical boundary- value 
solution near T,p in the form 

^ = - E E f dk °Mm ^t(t,r) ^j^r Y £m (Q) . (2.15) 
r fco m=-< - 7 - 00 sin^ij 

Here, the (nearly-) real quantities {akim} characterise the final data. 

An analogous description holds for fields of all the spins | , 1 and 2 that have so far 
been checked [2,5,7]. When considering perturbative boundary data for a field of any spin, 
posed on Ef in describing a final state resulting from black- hole evaporation, we denote 
by {a s k£mp} a set of analogous 'Fourier-like' coefficients, where s gives the particle spin, k 
the frequency, (£,m) the angular quantum numbers, and P = ±1 the parity (for s^O). 
For simplicity, we study in this Paper only bosonic perturbations, of spins s = 0, 1, 2, as 
treated in [2-4,6-8]. In each case, we found that the quantum amplitude or wave functional 
is of the semi-classical form, being given by 



{a s k£ m p}]T = N exp^iS^ {a sklmP } ; T ^ , (2.16) 



(2) 

where the pre-factor iV depends only on T . Here, S^iass denotes the (second- variation) 
action of the classical infilling solution, as a functional of the boundary data, corresponding 
to Eq.(1.2) for the s = case. For simplicity, we denote the collection of indices in 
a-sktmP by j . Further, we write Mj for the total (time-independent) ADM (Arnowitt- 
Deser-Misner) mass of the 'space-time', as measured at spatial infinity [27]. The ADM 
mass Mi , which is the limit at large radius of the variable mass m(t, r) of the Vaidya 
metric, is a functional of the final field configurations {cij} onEp, since it depends on the 
full gravitational field which results from finding the classical solution of the complexified 



8 



boundary- value problem. In Sec.6.1, we discuss this relationship between the total energy 
and the final field configurations, in the context of the normalisation of the probability 
density associated with the boundary data {aj} on Ep. 

As was found (for example) in the scalar case s = in [2-4,8], the classical action is 
dominated by contributions from frequencies k with the values 

TL 7T 

k = k n = — ; n = 1,2,3, ... . (2.17) 

It is also useful to define Akj to be the spacing between neighbouring /cj-values: 

Akj = 1 . (2.18) 

In an analogous way [2,7,25], for the corresponding Dirichlet problem with s = 1,2, 
and for s = (then neglecting the polarisation P ) , we found that 

o(2) 
class 



{a S k£mp} ; 

- h E t ± E * f^f *jf "*« < ~ wa M 



- ^M 7 T . 

(2.19) 

Equivalently, 

^class [{ a s^mp} ] T 

= ^2 X]( _1 ) SCs TiT~~\\ I dk k \ z sk£p\ 2 \a-sktmP + (~l) s P a s - k £ m p\ 2 cot(fcT) 
2 7 

(2.20) 

The parity operator P, taking the values ± 1 , is defined more explicitly through its action 
on the coefficients {a s k£ m p} , according to: 

a s k£mP = P (-l) m a>s,-ke,-mP ■ (2.21) 
The coefficients c s (s = , 1 , 2 ) are given by 

c = 2tt , ci = J , C2 = 1 . (2.22) 

For higher bosonic spins s = 1,2 in the adiabatic approximation above, the functions 
{£ s imp(t , r)} obey equations similar to the adiabatic version of Eqs. (2.4,5) for s = , but 
with a potential V s £p which depends on s : 

^ e (.b-a)/2 Q^j £ s£mp _ (d t ) 2 istmP + V siP ^simP = . (2.23) 



The explicit forms of V s £p for s = 1 and 2 are given in [2,7]. The complex coefficients 
{ z sktp} relate to the boundary conditions at spatial infinity for the radial part of the 
functions {^ s £mp(t,T)} , as in Eqs.(2.12,14) above for the case s = 0. Further details of 
the calculations for spins 1 and 2 are given in [2,7]. 

Eqs.(2.16,20) can be interpreted as giving a 'coordinate-representation' amplitude for 
each (square-integrable) final-field configuration specified by {a s k£ m p} > given that, on 
the initial hypersurface E/ {t = 0} , the perturbations vanish, or equivalently that the 
metric and background matter are spherically symmetric. Hence, \^[{a s k£ m p} ] T]\ 2 is 
the probability density for finding the field in a final configuration labelled by {a s k£ m p} 
at time separation T , measured at spatial infinity. 

(2) 

Following [2,7,25], the classical action functional S^iass * s f° un d to be a sum over 
individual 'harmonics' labelled by j , which depend on the corresponding indices {skjimP} 

2 

through the quantity \Aj\ , defined by 



2 (-1)' 



dj + (~1) S P ' O-s-kjlmP 



(2.24) 



The coefficients c s for bosonic spins s are given by Eq.(2.22). For s = 0, the quantities 
Zj are the complex numbers appearing in Eq.(2.14), which arise in solving the adiabatic 
radial mode equation (2.12); similarly for spins s = 1 and 2. This leads to the separated 
form of the quantum amplitude: 



{Aj} ; t] = N e -^ MlT l'[ (A,- ; T) 



(2.25) 



where N also depends only on T . 

Taking the classical action S^l ss in the form found in [2,25] for the scalar s = case 
(for example), one deduces that the wave functional for given scalar boundary data can be 
written as 



{Aj};T 



= N e" J 2 



\MjT -Q 



2isin(fc i T) 



exp 



-(A^)^|^| 2 cot(^T) 



This will be related to the coherent-state description in the following Section 3. 
3. Coherent states 



(2.26) 



The quantum amplitude (2.26) can usefully be re-written with the help of the Laguerre 
and Hermite polynomials [31], so as to demonstrate the relation with coherent states. First, 
we introduce the associated Laguerre polynomials Li m_n ' ) (x) , defined by 



r (m-n)( 



X) 



E 

p=0 



m \ (—x) 1 



n — p 



pi 



(3.1) 
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for m > n > . The usual Laguerre polynomials L n (x) [31] are given by 

L n {x) = £<?>(*) . 
The completeness relation for the set {L n (x)} reads: 



(3.2) 



J2 e" (x/2) L n (x) e"^ 2 ) L n (y) = S(x,y) . 



(3.3) 



n=0 



Writing z = x + iy , consider now the function L n (|z| 2 ) , which appears in the above 
quantum amplitude, re-expressed as in Eq.(3.5) below. For n > , this cannot be written 
as a product of two (decoupled) wave functions of x and y in an excited state, due to pair 
correlations [32]. But, in terms of Hermite polynomials H p (x) [31], one can expand 



L n (x 2 + y 2 ) = 



(~l) n 

2 2n n ! 



E 

p=0 



n 



P 



H 2p (x) H 2n - 2 p(y) ■ 



(3.4) 



With the help of the Laguerre and Hermite polynomials, the quantum amplitude 
(2.26) for our s = scalar boundary-value problem, arising in black-hole evaporation, can 
be re-written following Appendix A, in the form: 



{A,-};! 1 ] = Ne~^ MlT exp(-^ (Ak 3 ) k 3 \A 3 \ 2 /2) x 



x 



n e 



— 2iE n T r 



n=0 



kj (Akj) \Aj 



(3.5) 



where E n = (n + ^) kj is the quantum energy of the linear harmonic oscillator. Note also 
the dependence of the quantum amplitude on \Aj\ - it is spherically symmetric. 

There is a strong connection between the Schrodinger-picture wave functions 



N 



^ nj (xj,T) = - e~ ^/ 2) e" 2 ^" T L n ( Xj ) 



TV 



(3.6) 



appearing in the wave- function Eq.(3.5), where Xj = kj (Akj) \Aj\ 2 , and the exact quan- 
tum solutions to the forced harmonic oscillator [33]. In this approach, one considers a 
1-dimensional harmonic oscillator [33], with Hamiltonian 



H 



= h + ^"" v + gF(t) • 



(3.7) 



where F(t) denotes an external force, ji the oscillator mass and uj the oscillator frequency. 
Assume that F(t) =0 for t < t and for t > T , so that the asymptotic states, at early 
and late times t , are free-oscillator states. One can calculate the amplitude A^m to make 
a transition from the free-oscillator state \m > (with m particles) at early times t < to , 
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to the free-oscillator state \k > at late times t > T . Define the 'Fourier transform' of the 
force: 

(3 = [ dt F(t) e~ lLjt , (3.8) 

Jt 

and set 

* = (3.9) 
It has been shown [34-37] , in the case m > k , that 

- * a *- l ' m {0{^=T k *r>H*) , (3.10) 

where A is a real phase. This expression also gives A^m for m < k , since Ak m = A m k is 
symmetric. 

In the adiabatic limit, in which the force F(t) changes extremely slowly, one has 
z C 1 , and, from general considerations, a state which begins as \k > must end up in the 
same state \k > after the time-dependent force has been removed. From Eq.(3.10), one 
has 

A kk = e lX exp(-z/2) L k (z) . (3.11) 

The corresponding probability that there should be no change in the number of particles 
is |Afcfc| 2 = exp(— z) [Lk{z)] 2 . Apart from the introduction of mode labels j denoting 
the 'quantum numbers' {skimP} , together with a necessary re-interpretation for z , these 
amplitudes are effectively the wave functions Eq.(3.5) derived from our boundary- value 
problem. 

It will be useful to give a brief derivation of Eq.(3.10) in the context of the coherent- 
state representation. Coherent states \a > can be regarded as displaced vacuum states; 
that is, [13] 

\a> = D(a)\0> , (3.12) 

where 

D(a) = expect a) — a*a^j (3.13) 
is a unitary displacement operator, obeying 

D\a) = D-\a) = D{- a) , (3.14) 

and where the states \a > are eigenstates of the annihilation operator a with complex 
eigenvalue a . Among quantum states for the harmonic oscillator, they are the closest to 
classical states, in that they attain the minimum demanded by the uncertainty principle. 
Coherent states form an over-complete set, and are not orthogonal. In terms of the Fock- 
number eigenstates 

\n> = ^-L- |0> , (3.15) 
yn! 
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one has [13] 



\a > 



oo 

exp(-|a| 2 /2) Yl 



a' 



n=0 



\n > 



(3.16) 



The coherent state labelled by a = is the ground state of the oscillator. If, for example, 
the system started in a vacuum state, the amplitude to find it subsequently in a coherent 
state la > is 



< 0\a > - < 0|D(a)|0 > = exp(-|a| 2 /2) , 



(3.17) 



up to a phase. 

Proceeding towards a derivation of Eq.(3.10) by coherent-state methods, we note, 
from the properties of displacement operators, that 



£>(0I«> = £>(£)£>(«) |0> 

= exp[(V-r«)/2] D(£ + a)\0> 

= exp[(V-r«)/2] \£ + a> . 
For later reference, one can further show that 

D+( 7 )D^)D( 7 ) = D(») exp( 7 >- 7/ /) . 
Using Eqs.(3.16,18), one then has 



(3.18) 



(3.19) 



<m\D{£)\a> = 



(e + a) m exp 



-^(i«i 2 + i^i 2 + 2r«) 



(3.20) 



and 



a 



<m\D(0\a> = exp(-|a| 2 /2) ^ < m\D(0\n > 



n=0 



(3.21) 



On equating these, one finds that 



(1 + yr exp(-y|^| 2 ) = exp(|£| 2 /2) £ J^^y" <m\D{t)\n> 



n=0 



But, from the generating function for the associated Laguerre polynomials [31], 



(1 + yr e~ yx = 4 m ~ n) (*) 



\y\ < 1 , 



(3.22) 



(3.23) 



n=0 
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one deduces that the matrix element between initial and final states is 



< m\D(g)\n > 



ml 



l \ 2 



r~ n exp(-|^| 2 /2) 4 m "» )(|^| 2 ) , (3.24) 



which agrees with Eq.(3.10), up to an unimportant phase factor. 

4. Generalised coherent states 

The amplitudes < m\D(£)\n > , as in Eq.(3.24), can also be interpreted in terms of 
generalised coherent states |n, a > of the harmonic oscillator [35]. We define 



\n,a> = e 



■ iE n t 



D[a{t))\n> . 



(4.1) 



Then, in the Fock representation, one has 

oo 

\n,a> = < m\D^a(0) j \n > \m > 



iE m t 



(4.2) 



m=0 



As may be seen in Eq.(4.7) below, for generalised coherent states, the ground state (n = 0) 
is a coherent state and not a vacuum state. Generalised coherent states are to the coherent 
states what the Fock states \n > are to the vacuum state; that is, excited coherent states. 
In addition, one has, where I denotes the identity operator [35]: 



1 , V 

/ = — Ida \n,a> <n,a\ 

TV ' 



(4.3) 



<n,p\n,a> = L n (\a - /?| 2 ) exp^/Ta - \ (|a| 2 + |/3| 2 )^) , (4.4) 

d 2 a L n (\a - /3| 2 ) e r Q e ~^ 2 / 2 < n, o# > , (4.5) 



<n,/3|^> = 



e -|/3| 2 /2 



for an arbitrary state \ip > , with the definition: 



J d 2 a = J d Re (a) d Im (a) 



(4.6) 



Note that, from Eq.(4.4) with j3 = , one has 

<n,0|n,a> = <n\n,a> = e~' a ' ^ 2 L n (\a\ 2 j 



(4.7) 



again giving Eq.(3.11) up to a phase. The initial state could be seen not as a vacuum 
state, but as a Fock state, and the final state as a generalised coherent state. 

As shown by Hollenhorst [38], the amplitudes of Eq.(3.5) have yet a further interpreta- 
tion: they are the matrix elements for the transition from the state \k > to the state \m > 
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under the influence of a (linearised) gravitational wave, with the force F(t) proportional 
to the Riemann curvature tensor component R x txt(t) '■ 

F(t) = filR^it) = -\vt{dt?hH , (4.8) 

where t is the distance between two particles along the x-axis, each of mass ^ A* > while 
is the transverse-traceless gravitational-wave component of the metric [27] and x is 
the change in separation of the masses. 

In the context of black- hole evaporation, one expects that the role of the force is played 
by the time-dependent background space-time - which approximates a Vaidya space-time 
in the high-frequency limit at late times [2,6,20,21]. The above calculations indicate an 
explicit mathematical connection between the theory of forced harmonic oscillators and 
certain amplitudes relating to the dynamical evolution of black holes. 

An important point which we should mention is that, under the influence of a time- 
dependent force, an initial vacuum state transforms into a coherent state. Below, we 
discuss how, when one changes a phase parameter of the perturbations appearing in their 
frequencies (parametric amplification), an initial vacuum state transforms into a squeezed 
vacuum state. This phase is not an oscillator phase, but the small angle, 5 , through which 
the time T at infinity is rotated into the lower complex plane. 

5. Squeezed-state formalism 

In this Section and in the following Sec. 6, we shall see how, by rotating the asymptotic 
Lorentzian time T in the complex plane, and in the case of spherically-symmetric initial 
matter and gravitational fields, one obtains a quantum-mechanical highly-squeezed-state 
interpretation for the final state in black-hole evaporation, in the limit of an infinitesimal 
rotation angle. 

Grishchuk and Sidorov [16,17] were the first to formulate particle creation in strong 
gravitational fields explicitly in terms of squeezed states, although the formalism does 
appear in Parker's original paper on cosmological particle production [15]. In [16,17], it 
was shown that relic gravitons (as well as other perturbations), created from zero-point 
quantum fluctuations as the universe evolves, should now be in a strongly-squeezed state. 
Squeezing is just the quantum process corresponding to parametric amplification. 

Black-hole radiation in the squeezed-state representation was first discussed in [16,17]. 
The 'squeeze parameter' rj (see below) was there related to the frequency ujj and the 
black-hole mass M through 

tanh(rj) = exp(- 4tvM oo^j . (5.1) 

In this language, the vacuum quantum state in a black-hole space-time for each mode is 
a two-mode squeezed vacuum. However, our approach to squeezed states in black-hole 
evaporation is new, arising from a boundary-value problem involving two asymptotically- 
flat spacelike hypersurfaces, together with Feynman's +ie prescription [3]. We now give a 
brief account of quantum-mechanical squeezed states. 



15 



A general one-mode squeezed state (or squeezed coherent state) is defined [14] as 

\a,z> = £>(7)S(r,0)|O> = £>( 7 ) S(z) |0 > . (5.2) 
Here, .0(7) is the single-mode displacement operator, and we define 

S(r,</>) = S(z) = expQ(za 2 - z*a t2 )^ (5.3) 

in terms of annihilation and creation operators a and a* , respectively; we also define 

z = re~ 2i * . (5.4) 
Here, S(r,(f>) = S(z) gives the unitary squeezing operator for \a,z> , obeying 

S^z)S(z) = S(z)S^z) = 1 , (5.5) 

with 7 given by 

7 = a cosh r + a* e~ 2l< ^ sinh r . (5.6) 

The state Eq.(5.2) is a Gaussian wave-packet, displaced from the origin in position and 
momentum space. While the (real) squeezing parameter r ( < r < 00 ) determines 
the magnitude of the squeezing, the squeezing angle (\(f)\ < it/2) gives the distribution 
of the squeezing between conjugate variables. The squeezed vacuum state occurs when 
a = : 

\z> = \0,z> = S(z)\0> . (5.7) 

The limit of high squeezing occurs when r ^> 1 , for which the state \z > is highly localised 
in momentum space. 

Consider now the amplitude 

A = < a, z\D(n)\a, z > = < z\D^ '(7) D(fi) D(^)\z > . (5.8) 
One can use Eq.(3.19) to show that 

A = <z\D{n)\z> exp(^2zlm(7») , (5.9) 

and then use Eq.(5.6) and the relation a = \a\ e 1 ^ to show that [31] 

\A\ 2 = e" 17 ! 2 = exp^- |a| 2 (cosh2r + sinh 2r cos 2(0 + 0))^ . (5.10) 

Single-mode squeezed operators do not conserve momentum, since they describe the 
creation of particle pairs with momentum k . Two-mode squeezed operators, however, 
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describe the creation and annihilation of two particles (waves) with equal and opposite 
momenta. A two-mode squeeze operator has the form [38] 



S(r,</>) = expl r[ e~ 2 ^ a+ a_ - e 2t(p a\ al ) I , (5.11) 




where a± and a± are annihilation and creation operators for the two modes, respectively. 

Consider two conjugate operators, p and q , with variances Ap and Aq . In the 
squeezed- state formalism, one may construct states such that Ap and Ag are equal, taking 
the minimum value possible. The name 'squeezed' refers to the fact that the variance of 
one variable in a conjugate pair can go below the minimum allowed by the uncertainty 
principle (the squeezed variable), while the variance of the conjugate variable can exceed 
the minimum value allowed (the superfluctuant variable) [14,40,41]. The superfluctuant 
variable is amplified by the squeezing process, and it becomes possible to observe this 
variable macroscopically; in contrast, the subfluctuant variable is squeezed and becomes 
unobservable. In particle production, whether by black holes or in cosmology, the number 
operator is a superfluctuant variable, while the phase is squeezed. 

6. Analytic continuation and the large-squeezing limit 

We return to the quantum state described in the Schrodinger picture by Eq.(2.20), 
giving the wave-functional \I/[{Aj}; T] for perturbative bosonic field configurations on 
the final surface T, F , labelled by 'coordinates' Aj = A s g m p . Again, T denotes the 
(Lorentzian) time separation, measured at spatial infinity, between the initial surface Ej 
and the final surface Hp ■ For convenience, we repeat Eq.(2.26) and Eq.(3.5): 



= Ne~^ MlT J] y(Aj;T) 
j 

s * <-«"•* n {Aki) k > lAif cot{k ' T) ) 

oo 

= Ne-* M > T exp(-^ (Akj) kj |A,| 2 / 2 ) ]J ^ e~^ T L n ( kj (Akj) |A,| 2 ) . 

j j n=0 

(6.1) 

We now define 

{A,}; T] = Ne~^ M ' T \[ ^sin^T) ^(A,; T) 

j 

ee N e~^ M ' T exp(l (Akj) kj |AJ 2 cot(^T)] ( 6 .2) 

j ' 

= iVexp^2 ss [{A,};T 
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where N is a T-dependent prefactor. 

As described in the Introduction, the classical boundary-value problem is expected 
to become well-posed when one rotates the asymptotic time-interval T into the complex, 
taking 

T = \T\e~ iS ; < 5 < - . (6.3) 



By contrast, if one did not rotate T , so that T remained real and positive, then the 'sum' 
in Eq.(6.2) would diverge, due to the simple poles on the real-frequency axis at 



n = 1,2, ... 



(6.4) 



assuming that kj \Aj\ remains suitably non-zero near kj = a n . Thus, the quantum am- 
plitude cannot be computed simply by working with space-times of Lorentzian signature. 
As in the Introduction, we follow Feynman [1] in adopting a + is prescription. That is, we 
carry out the quantum calculation above for 5 > , and then, at the end of the calculation, 
take the limit 5 — > + of the quantum amplitude, to obtain the Lorentzian amplitude. 

Note also that we will not eventually take the limit of infinite |T|. However, we do 
expect that, in practice, |T| will far exceed the dynamical collapse time-scale for the black 
hole, which is of order it Mi [27]. We shall see that, in our problem, the replacement 
Eq.(6.3) of real T by T = \T\ exp(— id), < 5 < tv/2 , leads to a squeezed state, with 
a high degree of squeezing for small 5 . For a general 5 > , one has, from Eqs. (6.2,3): 



{Aj};T 
= $ 



x 



{Aj}; \T\,S 



-i|Mj|T|cos(5 - iM 7 |T|sin(5 



X 



| J expf — - (Akj) kj \Aj\ 2 coth^/cj \T\ sind — i(pj(\T\ , 5)j j 
j ' 

N g-IO^TIcos* e -*M,|T|sin* "Q 6 X P ^ (of } + Z Of } ) (Afc,) ^ | ^ 



Here, we define 



>f(l 



— kj |T| cos 5 , 



ni R) ( \t\ , <s) 



sinh^/cj |T| sin 5) 



nfi\r\,<) 



2 ^cosh 2 (/cj \T\ sin 5) — cos 2 4>j^j 
sin (2 (f)j) 
2 (cosh 2 (kj \T\ sin 5) — cos 2 0jj 



(6.5) 

(6.6) 
(6.7) 

(6.8) 
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One can further re- write Eq.(6.5) in the form: 



d>[{^}; \T\,8_ 

- N e - ^Mj^coss - ±Mj\T\sin6 TT exT Y_ ]_ ( Ak .) / 1 + e 2 ^ tanhr A , 
" 7V 6 6 11 6XP V 2 [/Xk i> kj \l- e^tanhrj ' ■ 



,2 



where we have set 



(6.9) 

tanh n (\T\,6^ = exp(-2/c i \T\ sm5) . (6.10) 



Therefore, 



exp(-2r i ) = tanh(fcj \T\ sintf) . (6.11) 



Hence, rj — > for high frequencies, while — > oo for low frequencies. On comparing 
with Sec. 5, we recognise Eq.(6.9) as the coordinate-space representation of a quantum- 
mechanical squeezed state [23,42], with r,-(|T| , 5) the squeeze parameter and 4>j(\T\ , 6) 
the squeeze angle. The evolution of the squeezed state is taken into account by the 
\T\ dependence in rj and in <pj , which are in general complicated functions of time. 
Eq.(6.6) becomes simpler in the limit of infinitesimal 5 . Neglecting terms of 0(5 2 ) , one 
has 4>j{\T\ , S) ~ — kj \T\ , corresponding to free evolution. 

Computing the probability density ; \T\ , S]\ , one finds that 



*[{^-}; \t\,s 2 



cothe, , . . x . , . , 2 \ ( 6 - 12 ) 



= |iV| 2 exp(-M 7 |T| sin s) ]J exp (- /(fc ^ ^ ) (Afc 3 -) fc 3 - |A 
where 

/fc,6,-,|T|) = 1 + Sm fa |r|) , (6.13) 
V / sinh e 3 

and 

ej = fcj|T|sin<J . (6.14) 

Eq.(6.12) describes a Gaussian non-stationary process, in that the variance is an 
oscillatory function of time. One is now dealing with standing bosonic waves, rather than 
with travelling waves; the (classical) amplitudes for left- and right-moving waves are large 
and almost equal, much as in the cosmological scenario [16,17]. These standing waves 
imply a correlation between particles with opposite frequencies (and azimuthal angular 
momenta ± m ) in the final state. One consequence of the high-squeezing behaviour is that 
the variance for the amplitudes {Aj} is large, so that there are large statistical deviations 
of the observable power spectrum from its expected value. This is just a manifestation of 
the Uncertainty Principle. 
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We now assume that 5 is sufficiently small that tj <C 1 , or, equivalently, that 
< S < (kj |T|) _1 . Then, from Eq.(6.11): 



€j ~ exp(— 2rj) 



ej « 1 



(6.15) 



corresponding to r, ^> 1 , which is the high-squeezing limit. Thus, the high-squeezing 



limit corresponds to the limit 5^0+ with kj 



\T\ bounded 
.(«) 



'] 



A broadening of the width 
, therefore corresponds to a large 



of the position distribution, which is of order [fi 
squeezing in the momentum distribution. 

In the squeezed-state formalism, we regard this as the classical limit, since the average 
number of particles in the final state is large: 



<Nj> 



= 4 sinh 2 



exp (2 Tj) 



(6.16) 



for rj 3> 1 . For another way in which to view this, consider again the state Eq.(6.5). 
The WKB condition is met when 



sin(2 0j) 



sinh(2ife j \T\ sin S) 



> 1 



(6.17) 



which is satisfied in the high-squeezing limit e 3 <^ 1 . The final state of the remnant black- 
hole evaporation flux, therefore, becomes more classical in the WKB sense in the limit 
5 — > . In this limit, one can effectively consider the final perturbations as being repre- 
sented by a classical probability distribution function [16,17,40,43]. As in the inflationary 
scenario in cosmology, the perturbations away from the spherically-symmetric black-hole 
background space-time, of quantum-mechanical origin, cannot be distinguished from clas- 
sical stochastic perturbations, without the need of an environment for decoherence. 

An inflationary analogue can also be described for the initial conditions on the per- 
turbations in the black-hole case. In cosmology, the assumption is that, at some early 
time just prior to inflation, the modes are in their adiabatic ground state. This, in turn, 
originates from the assumption that the universe was in a maximally-symmetric state at 
some time in the past [18]. A similar assumption is present in our black-hole case, where 
we assumed that the initial perturbations were very weak, so that the initial 'star' and its 
gravitational field were spherically symmetric. 

To obtain one further view on the late-time state, for small ej , one can express 
Eq.(6.12) in the form 



*[{^-};|T|,<5->0; 



|Af axp(-M,|rlJ)n«xp(- (ei)a+ ^ fe lrl) 
|iV| 2 exp(-M/ \T\ 5) Yl exp^- (Akj) kj \A 3 \ 2 p^j , 



(Akj) kj \Aj 



(6.18) 
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where we have defined 



oo oo 

Pj = (Aw n ) 8(kj - oj n ) = ^ exp(2mfcj \T\) , (6.19) 

n= — oo n= — oo 

and used the delta-function identities 

*(*) = i hm ' , (6.20) 

7T e^O [e z + X Z ) 



and 



'(/(*» = ^TFRT ' (6 ' 21) 



where the Xi are defined to be the zeros of f(x) . In these equations, we are also using the 
definitions u n = mr/\T\ and Auj n = iv n+ i — u n = ir/\T\ . Hence, interchanging the 
sums over j and n , in the continuum limit for the {kj} frequencies, one has 



2 'tmax 

^[{Aj}; |T|,<?-0]| ~ |iV| 2 J] J] exp(-(Aa; n )a; n |^ m p| 2 ) , (6.22) 

simP n=l 



in the small-5 limit. The sum Xl^L-oo nas been converted into J2^=i > noting that 
fcj > 0, and also that kj \Aj\ — >0 as kj — »0, and we have introduced n max , the largest 
value of n such that cu nmax = A^j , providing an effective cut-off in the product over n . 
This agrees with the result summarized in Eq.(19) of [3], where the calculation of the 
probability density was based on the contour-integration treatment of [2]. 

The presence of the delta function in Eq.(6.18) indicates that, in the high-squeezing 
limit, the random variable <f>j associated with the final state is squeezed to discrete values, 
independently of the quantum numbers {s£mP} . Note that it is only the squeeze phases 
{4>j} of the (standing- wave) perturbations which are fixed and correlated in the high- 
squeezing limit. 

6.1 Normalisation 

We now discuss the normalisation factor \N\ in the probability density Eq.(6.12). 
We consider the dimensionless variables {xj} , as defined after Eq.(3.6). Now, |iV| 2 is 
determined by integrating the probability density Eq.(6.12) over the space of all {xj} , 
since the sum of all probabilities of all possible configurations {xj} is unity. Hence, 



,, T ,2 -r-r coshe, sinhe 

\N\ =11 



(sinh 2 ^ + sm 2 (kj |T|)) 

n 



(6.23) 



cosh(2rj) — cos (20 j) sinh(2r 
One can verify from Eqs.(6.6,ll) that this infinite product converges 
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There is, however, an ambiguity in the form of the normalisation factor due to the 
presence of the surface term, — | Mj T , in the action Eq.(2.20), arising from the boundary 
at large radius joining the initial and final space-like hypersurfaces. The origin of this 
ambiguity lies in the fact that the total ADM mass Mi is a functional of the final field 
configurations {xj} - see Eq.(8.6) below. 

Thus, in this case, 



\N\ 2 



J] (A 2 Of) +26,) , (6.24) 



where we have re-introduced the factor of A 2 : here, A is the small parameter used in the 
expansion of the metric and fields in Eq.(l.l). Ambiguities caused by surface terms in the 
action were also discussed in the squeezed-state formalism in [42]. 

One consequence of the high-squeezing behaviour is that the variance for the ampli- 
tudes {Aj} is large, so that there are large statistical deviations of the observable power 
spectrum from its expected value. This is just a manifestation of the uncertainty principle. 
Indeed, from Eq.(6.23) and the probability distribution Eq.(6.12), we find that the average 
value of Xj has the form 



< x j > 



= (cosh(2rj) - cos(20 i )sinh(2r J -)) . (6.25) 

Thus, in the large-squeezing limit r j ^> 1 , we have 

< Xj > ~ ^ exp(2r-,-) sin 2 (<pj) . (6.26) 
From Eqs.(6.14) and (6.15), one finds 

1 sin 2 (6 A 

<x ' > = y W ' (6 ' 27) 

It is clear from this equation that < Xj > is not a smooth function of j , since it has 
a large number of peaks, and also has zeros when <pj = nn , for each integer n . This is 
indicative of the standing-wave feature of our boundary-value problem. 

For comparison, in inflationary cosmology the oscillation phases of standing waves 
have fixed values, giving rise to zeros in the power spectrum, which are characteristic of 
the CMBR. The power spectrum of cosmological perturbations in the present universe is 
not a smooth function of frequency. The standing- wave pattern, due to squeezing, induces 
oscillations in the power spectrum. This in turn produces Sakharov oscillations [43,44], 
due to metric and scalar perturbations, in the distribution of higher-order multipoles of 
the angular correlation function for the temperature anisotropies [24,45] in the CMBR, 
for all perturbations at a given time whose wavelength is comparable with or greater than 
the Hubble radius defined for that time. That is, the peaks and troughs of the angular 
power spectrum have a close relationship with the maxima and minima of the metric power 
spectrum. For long wavelengths, the power spectrum becomes smoother. 
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7. Entropy and Squeezing 

There have been many accounts of how to determine entropy generation in the squeez- 
ing formalism [32,39,42,43,46-49]. Hu and Pavon [50] were the first to associate entropy 
generation with the monotonic increase in particle number with time, induced by para- 
metric amplification in a vacuum cosmological space-time. As squeezing is the quantum 
analogue of parametric amplification, one would expect that entropy production could be 
calculated via the squeezed-state formalism. This is indeed the case, although, as with any 
entropy calculation, the nature of the coarse-graining must be specified. For squeezing, 
this is particularly relevant since squeezed evolution is unitary: that is, there is in principle 
no loss of information in the evolution of the initial pure state to the final pure squeezed 
quantum state. The definition of entropy depends on how one chooses to measure the ob- 
servables associated with the final squeezed states. For example, one can reduce the final 
density matrix with respect to a Fock or coherent-state basis [42], or use eigenstates of 
the superfluctuant variable. In [32,40,41], the loss of information comes from the increased 
dispersion of the superfluctuant operator. 

Following the work of [32,46-49], a universal form for the entropy density growth Ar- 
nolds for each mode, when one studies the classical limit of large average particle number, 
corresponding to the large-squeezing regime, namely 

ASj ~ 2r 3 ; Tj > 1 , (7.1) 

irrespective of the particular coarse-graining. On calculating the von- Neumann entropy S 
from Eq.(6.12), using also Eq.(6.23), one finds that (in units such that ks = 1): 

S = - J \[dx j P{ Xj ) \ogP{ Xj ) 

= 1+2-/ logfe 2r ''sin 2 0j- +e 2rj cos 2 <j>A . 
j 

Thus, the entropy (7.2) arises from our ignorance of the precise final radiation configura- 
tion. In the high-squeezing limit, one has, therefore, 

AS ~ 2^ Tj + J2 lQ g( sin2 ^) • ( 7 - 3 ) 

j j 

which agrees with Eq.(7.1) when rj 3> 1, even if sin^j ~ . Eq.(7.3) remains valid 
even if we take into account the ambiguity in the normalisation factor discussed in Sec.6.1. 

8. Classical predictions 

We now discuss the way in which strong peaks in the wave function lead to some 
definite predictions. In quantum cosmology, wave functions ^(q) are commonly peaked in 
such a way as to describe (semi-classically) families of classical trajectories, corresponding 
to solutions of the Hamilton-Jacobi equation. Loosely speaking, such wave functions are 
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peaked about correlations between coordinates and momenta. Such correlations may per- 
haps be described more clearly with the help of the Wigner function W(p, q) [51]; in this 
reference, the correlations were precisely identified via W(p, q) . The opposite extreme - 
no such correlation - occurs when W(p, q) factorises into a product of a function of position 
q and a function of momentum p . On the other hand, when W(p, q) is peaked about some 
'surface' in phase space, say {p = /(<?)}, then the wave function predicts this particular 
correlation. 

An alternative proposal for measuring correlations was given in [52]. There, classi- 
cal correlations for a given Wigner function were predicted by means of projection onto 
coherent states, where position and momentum are equally (un) known, as in the classical 
theory. As demonstrated in [51], in the harmonic-oscillator case, the correlation between 
p and q is such that the Hamiltonian equals the classical energy. In the present paper, we 
arrive at similar conclusions. 

We now return to our wave function #[{A,};T] of Eqs.(2.25,26), describing the 
quantum amplitude for typical anisotropic perturbations of spins s = , 1 and 2 about 
a final background spherically-symmetric gravitational and massless-scalar field (7^ , $) . 
This amplitude is further described in Eq.(3.5), in terms of a product over modes j , each 
term involving an exponential times a Laguerre polynomial. The relevant argument of each 
Laguerre polynomial is the dimensionless quantity xj , given after Eq.(3.6). Following the 
above discussion in this Section, we look for predictions from the Heisenberg-picture wave- 
functional, given by a product over j of the wave functions 

= ^ e^(-Xj/2)L nj ( Xj ) . (8.1) 

If we were to restore all units, the dimensionless argument of the Laguerre polynomial 
would acquire a factor h~ x . In most cases of astrophysical interest, one then has x j 1 . 
Examining the wave function (8.1) in the limit of large argument, note from [31] that 

(— x) n 

L n (x) ~ : — , as x — > 00 . (8.2) 

n! 

For each j , one can find the peak in the wave function as a function of Xj , which is at 

1 



Tin = 



kj (Akj) \Atf h- 1 . (8.3) 



3 2 

Taking the spin-0 case, for example, this gives, from Eqs. (2.24,8.3), 



n 



3 



2tt kj (Akj) \zj\ aj + ao-k^mp . (8.4) 



2 



But in Sec. 2 of [8] (see also [4]), we showed in Eq.(2.17) that, for spin-0 perturbations, one 
has 

(8.5) 



2 9 ^ 

\bj\ = 2n kj \Zj\ 



ttj + a kilmP 
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where the {bj} are 'Fourier amplitudes' associated with the radiation reaching X + (future 
null infinity). Thus, Eq.(8.5) gives a match between the positive- frequency decomposition 
for particles reaching X + (travelling waves) and the boundary- value formulation employed 
in our papers [4,8], with final field configurations specified on the hypersurface T,p , given 
by {t = T} (standing waves). From Eqs.(8.4,5), one has nj = (Akj) \bj\ , for each j . 
Thus, we find that 

Y^hnjkj = (Akj) kj \bi\ 2 = M I . (8.6) 

j j 

The left-hand side is just the total energy of the radiated particles. The middle expres- 
sion is the total energy in the massless-scalar fluctuations, which (in the absence of any 
gravitational radiation) equals the initial ADM mass Mj . 

9. Conclusion 

In this paper, we have illustrated many aspects of the boundary-value formulation 
for linearised integer-spin fields propagating in an evaporating black-hole space-time. For 
simplicity, we have taken here only the case in which one has initial fields which are 
spherically symmetric. When the (Lorentzian) proper-time separation T at spatial infinity 
between the initial and final hypersurfaces is deformed infinitesimally into the lower-half 
complex plane, following Feynman [1], one obtains a quantum-mechanical squeezed-state 
formalism. The large-squeezing limit is equivalent to the WKB limit, and corresponds to 
rotating T by only an infinitesimal angle 5 into the lower half-plane. Since the final highly- 
squeezed state is a pure state, the unpredictability associated with any final momentarily- 
naked singularity in the Lorentzian space-time is avoided. 

We found that, as in the cosmo logical case, the bosonic perturbations on the black- 
hole background can be regarded as a stochastic collection of standing waves, rather than 
as travelling waves, in the high-squeezing limit. This leads to the prediction of peaks in 
the power spectrum, for 'relic' radiation from an evaporating black hole, analogous to the 
Sakharov oscillations in the CMBR. 
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Appendix A: Derivation of the wave function Eq.(3.5) 

Here, we describe the chain of relations leading, in the scalar case, from Eqs. (2.25,26) 
to Eq.(3.5), thus expressing the quantum amplitude \I/[{A,}; T] as an infinite product, 
over the index j , of suitable exponentials and Laguerre polynomials. We first set 
a = exp(— ikjT) , y = and x = Xj in Mehler's formula [31] 

(!-„»)* f; e -*CW) W)g>(») = e i<*'-»'>exp(-(*- 0! ,)V(l-a 2 )), 



2Pp\ 



(A.l) 



where 



H p (x) = (-If e* 2 ^ (e-* 2 ) ; p = , 1 , 2 , . . . (A2) 
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define the Hermite polynomials, which appeared already in Eq.(3.4). This normalisation 
of the H p (x) implies that 



2Ppln 

With the above choice of a , y and x , we find 

1 



r J dx exp(-x 2 ) \H p (x)j = 1 . (A3) 



(2m sin(fcjT)) exp [ - i (xj) 2 cot (kjT) I = ^ exp^— zE^T^ ip p (xj) ip p (0) , 

(A4) 

where E p = kj [p + |) , and we define 



exp(- 2 (^) 2 ) ifpfo, 
Vv(^) = J — , (A.5) 

such that 

/oo 
cfc,- |V P (^)| 2 = 1 • (A.6) 
-oo 

Only even terms contribute to the sum in Eq.(A.4), since [31] 

H , 2p{0) = turn , ( ,. 7) 

pi 

H 2p +i(0) = . (A8) 

Hence, 

{2-ni sm(kjT) ) expf -i (xj) 2 cot (kjT) J = exp(- i E 2p i/j 2p (xj) ifj 2p (0) . 

(A9) 

In addition, we now define, for each index j in Eq.(3.5), suitably scaled versions of the 
real and imaginary parts of the complex quantity Aj by: 



x 3 = [ k A 2 Re(Aj) , Vj = (&\ 2 i m{Aj ) . (A.10) 
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Then, 

exp^ikj V~ x \ Aj\ 2 cot(fcjT) 
2m sm(kjT) 

oo oo 

= ^2 ex p(~^(^2p + E 2p ')T\ lp2p{Xj) fop'iVj) ^2p(0) ^2p'(0) 

p'-0 p=0 
oo p 

= ^ GXp(-2iE p T) ^ ^2p' (Xj) 1p2p-2p' (Vj) VV(°) ^2p-2p' (0) 
p=0 p'=0 



7T 



1 exp(-fcj |Aj| 2 /2y) x 

~ (_i)P exp(-2»E p r) ^ gl 

2>^! ^ (p')\(p-p')\ H2p ' {Xj) H *»--#KVi, 



p=0 p'=0 

(All) 

where we define 1/ 1 = Akj . In the limit \T\ — > 0_|_ (or rather, fcj |T| <C 1 ), since the 
eigenfunctions ip p (xj) form a complete orthonormal set, one has 

lim (2ni sin (fcj T) j exp^zfcj |A,-| 2 cot(fcjT)/2T^J = S(xj) S(yj) 



(A12) 

This suggests that Aj — ► as fcj \T\ — > 0+ , agreeing with our initial condition that the 
fields are spherically-symmetric. We now use the identity Eq.(3.4) in Eq.(A.ll). Taking 
the product over all j , and introducing a normalisation factor N , we find 

*[{^-};t] = n^(^; T ) 

j 

= Ne-^ T exp(-5> \A//2V) ]J f] e" 2 ^L p (^ |A,f) , 

(A13) 

which gives Eq.(3.5), where we have now included the contribution to \1/ from the time-like 
boundary near spatial infinity One can confirm that Eq.(A.13) also gives Eq.(2.26); this 
simply involves use of the generating function for Laguerre polynomials [31]. 
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